data=read.table('bac-2008-02-13.blast')
pdf("bacplot10.pdf")
hist(data$V3,breaks=10)
hist(data$V4,breaks=10)
hist(data$V5,breaks=10)
hist(data$V6,breaks=10)
hist(data$V7,breaks=10)
hist(data$V8,breaks=10)
hist(data$V9,breaks=10)
hist(data$V10,breaks=10)
hist(data$V11,breaks=10)
hist(data$V12,breaks=10)
plot(data$V12,data$V11)
plot(data$V4,data$V11)
plot(data$V4,data$V12)
q(save="no")

